Unconventional metamagnetic electron states in orbital band systems 
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We extend the study of the Fermi surface instabihty of the Pomeranchuk type into systems with 
orbital band structures, which are common features in transition metal oxides. Band hybridiza- 
tion significantly shifts the spectra weight of the Landau interactions from the conventional s-wave 
channel to unconventional non-s-wave channels, which results in anisotropic (nematic) Fermi sur- 
face distortions even with ordinary interactions in solids. The Ginzburg-Landau free energy is 
constructed by coupling the charge-nematic, spin-nematic and ferromagnetic order parameters to- 
gether, which shows that nematic electron states can be induced by metamagnetism. The connection 
between this mechanism to the anisotropic metamagnetc states observed in Sr3Ru207 at high mag- 
netic fields is studied in a multi-band Hubbard model with the hybridized quasi-one dimensional 
dxz and dyz-bands. 

PACS numbers: 71.10.Ay, 75.10.-b,75.30Kz 



I. INTRODUCTION 



Pomeranchuk instabilities are a large class of Fermi 
surface instabilities in both density and spin channels, 
each of which can be further decomposed into differ- 
ent partial wave channels^. This class of instabilities 
in the non-s-wave density channel result in uniform but 
anisotropic (nematic) electron liquid states^, which have 
recently attracted a great deal of attention in recent 

■■.^,^j. j2. 3 .4. 5. 6. 7. 8. 9 . 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21 ^xi partic- 
ular, these instabilities have been studied in the con- 
text of doped Mott insulators^^, high Tc material a'^'^^ , 
and quantum Hall systems in nearly half-filled Landau 
levels^'^'^'*. Experimental evidence has also been found 
in ultra-high mobility two-dimensional electron gases and 
quantum wells in large magnetic field o^^'^^i^^ . 



Non-s-wave spin channel Pomeranchuk instabilities 
are ^^unconventional magnetism" in analogy to uncon- 
ventional superconductivity, which have been exten- 
sively investigatedMi5i28^^30^,32^ j^^ ^^^^ by ^u 

and Zhang, these states are classified as isotropic and 
anisotropic phases dubbed (3 and a-phases as the coun- 
terparts of ^Yie-B (isotropic) and A (anisotropic) phases, 
respectively. The /3-phases have circular or spherical 
Fermi surfaces with topologically non-trivial spin con- 
figurations in momentum spaced. In the /3-phase, the 
relative spin-orbit symmetry is broken, a concept intro- 
duced in the context of ^He-B phase, while essentially 
the overall rotational symmetry is not. The a-phases are 
anisotropic electron liquid crystal states with spin degree 
of freedom, which have been studied by many groups: 
the p-wave phase was first studied by Hirsch^>22 under 
the name of the "spin-split" state, and was also proposed 
by Varma et al^^ as a candidate for the hidden order 
phenomenon in the heavy fermion compound URu2Si2; 
the d-wave phase was studied by Oganesyan et al.^ under 
the name of "nematic-spin-nematic" phase. Systematic 
studies of the ground state properties and collective exci- 
tations in both the anisotropic a and isotropic /3-phases 
has been performed by Wu et al^^^. Very recently. 



Chubukov and Maslov found that when approaching the 
ferromagnetic quantum critical point, the p-wave chan- 
nel spin Pomeranchuk instability develops prior to the 
developing of ferromagnetic instability^. 

Although unconventional magnetism has not been con- 
vincingly identified in experiments, a spontaneous ne- 
matic electron liquid has been observed in the ultra-clean 
samples of the bilayer ruthenate Sr3Ru20 , ■^"^i^^iSe ^ which 
has arouse much research interesl4 '^^i'^^'^^i'*°''^^i"^^'^^ . 
Sr3Ru207 is a metallic itinerant system with the tetrag- 
onal Ru02 (ab) planes. It is paramagnetic at zero mag- 
netic field, and develops two consecutive metamagnetic 
transitions in the external magnetic field B perpendic- 
ular to the a6-plane at 7.8 and 8.1 Tesla below IK, re- 
spectively. In the state between two metamagnetic tran- 
sitions, the resistivity measurements show a strong spon- 
taneous in-plane anisotropy (nematic) along the a and b- 
axis with no noticeable lattice distortions, thus this effect 
is of electronic origin. It can be interpreted as a nematic 
state with the anisotropic distortion of the Fermi surface 
of the majority spin polarized by the external magnetic 
field in Rcf.'^^, which is essentially a mixture of the d- 
wave Pomeranchuk instabilities in both density and spin 
channels. 

In spite of years of intensive research, most theories 
remain phenomenological without considering the ac- 
tual orbital band structures of Sr3Ru2 07^i^i^. Two 
key questions have not been answered satisfactorily after 
years of intensive research. First, Sr2Ru307 is a t2g-heLnd 
system, having both quasi-one dimensional bands of d^z 
and dyz and the two-dimensional bands of d^y. Which 
are responsible for the nematic behavior? Second, the ne- 
matic states require strong exchange interactions in the 
d-wave channel, but the usual exchange interaction from 
Coulomb repulsion is mostly in the s-wave channel. Mi- 
croscopic theories based on the single band picture of dxy 
combined with the van Hove singularity have been nicely 
developed^i^i^. However, their models need a large d- 
wave channel exchange interaction which was introduced 
by hand. It is difficult to justify the microscopic origin 
of this interaction in terms of Coulomb interaction. Fur- 
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thermore, these theories do not address the fact that the 
nematic ordering does not appear in the single layer com- 
pound Sr2Ru04, which has a similar band structure of 

In this article, we provide a natural answer to these 
two key questions by extending the theory of Pomer- 
anchuk instabilities into multi-orbital systems. We point 
out that it is the hybridized quasi-one-dimensional d^z 
and dyz bands instead of the d^y band that are responsi- 
ble for the nematic ordering based on the following rea- 
soning. The key difference of electronic structures be- 
tween Sr3Ru2 07 and Sr2Ru04 is the bilayer splitting, 
which is prominent for the quasi-one dimensional bands 
of dxz and dyz but small for the two-dimensional bands 
of dxy It is natural to expect that the spontaneous ne- 
matic behavior occurs in the bands of d^z and dyz and 
is accompanied by an orbital ordering. Furthermore, the 
orbital band hybridization between them shifts a signifi- 
cant spectra weight of the exchange interaction into the 
d-wave channel, thus the nematic ordering can arise from 
the conventional multi-band Hubbard interactions. This 
mechanism also applies to other strongly correlated or- 
bital systems. 

This paper is organized as follows. We first present a 
heuristic picture to illustrate the idea how orbital hy- 
bridization enhances the Landau interaction in non-s- 
wave channels in Sect. |TT1 In Sect. IIIIl a phenomcnolog- 
ical Ginzburg-Landau free energy is constructed to ex- 
plain the two consecutive metamagnetic transitions and 
the nematic phase in between. In Sect. IIV| we use a 
microscopic multi-orbital Hubbard model based on the 
quasi-one dimensional dxz and dy^-bands to explain the 
appearance of the nematic state. Conclusions and out- 
looks are made in Sect. |Vl 

II. LANDAU INTERACTIONS MODIFIED BY 
ORBITAL HYBRIDIZATION 

In this section, we present a heuristic picture to illus- 
trate the enhancement of the non-s-wave Landau interac- 
tions from orbital band hybridizations. For a single band 
system without orbital structures, the Landau interac- 
tion functions can be simply expressed at the Hartree- 
Fock level in the density and spin channels as: 

r(pi,P2) = V{q^0)^^V{\pi-p2\), 

riPuA) - -impi-P2i), (1) 

where V{p) is the Fourier transform of the two-body in- 
teraction T^d^i —'^2 1), say, the Coulomb interaction. The 
high partial wave channel components of V are usually 
weak, thus the condition for Pomeranchuk instabilities in 
high partial channels is more stringent than that of the 
s-wave instability of fcrromagnetism. 

This situation is significantly changed in multi-band 
systems with non-trivial orbital hybridization. Let us 



consider a simplified two dimensional example of the hy- 
bridized bands between dxz and dyz and assume the sin- 
gle particle eigenstates around the Fermi surface takes 
the form of 

\Mp}) = e^nuip))®X<r, 
\u{p)) = COS (f>p\dxz) + sin (fipldyz) (2) 

where (j)p is the azimuthal angle of p; \u{p}) is the Bloch 
wave function with the internal orbital configurations; 
Xcr((T =t, I) are the spin eigenstates. This orbital struc- 
ture has no effects on the Hartree interaction between 
two electrons with opposite spins, while it significantly 
changes the Fock exchange interaction of two electrons of 
the same spin which are not completely indistinguishable 
any more. Consequentially, the exchange interaction be- 
tween them acquires an extra form factor describing the 
inner product of their orbital configurations as 

-m-P2)iwpi)KP2))p. (3) 

The Landau interaction functions change to 

riPuP2) = 1/(-f-O)-i[l + cos2(0p, -0pj] 

X V{pi -P2), 
nPi,P2) = -i[l + cos2(0p, -0p,)]y(pi-p2). (4) 

Therefore even if the bare interaction V{pi — P2) is 
dominated by the s-wave channel, the extra d-wave factor 
arising from the orbital hybridization shifts a significant 
weight into the d-wave channel. 

III. THE GINZBERG-LANDAU THEORY 

The Eq. |4] above implies that a strong ferromag- 
netic or metamagnetic (s-wave) tendency also enhances 
the nematic ordering (d-wave). Before constructing the 
microscopic theory, we illustrate this point through a 
Ginzburg-Landau free energy formalism which includes 
the coupling among the ferromagnetic order m, the 
charge nematic order Uc, and spin nematic order Usp. 

In the square lattice, the two different nematic chan- 
nels dx^~y2 and dxy belong to non-equivalent represen- 
tations. Only the dx^-y^ channel instability is experi- 
mentally observed and thus is kept. Due to the experi- 
mentally observed anisotropy between the z-axis and the 
a6-plane, we only keep the z-component of spin and spin- 
nematic orders. The Ginzburg-Landau free energy is con- 
structed as 

F{h) = F{m) - hm + rcul + rspulp + gcul + gsn*p 
+ g{m)ncnsp (5) 

where F{m) is the magnetic order contribution to the 
free energy as an even function of m; h is the external 
magnetic field; Tc oc (2 -t- F2) and r^p oc (2 -I- F2) are the 
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FIG. 1: a) The sketch of the magnetic part F(m) in the 
Ginzburg-Landau free energy. The slopes of two common 
tangent hnes mark the magnetic fields of meta-magnetic tran- 
sitions, b) The dimensionless coupling function g'irn) has 
a peak distribution between two consecutive matamagnetic 
transitions. 

mass terms of charge and spin nematic orders, respec- 
tively; nc,sp are defined as: 

k.a 

risp = ^{■ilji{k)a^.ai3tppik))f{k), (6) 

where f{k) is a form factor exhibiting the dx2_y-2 sym- 
metry; g{m) is the coupling function between Uc and Ugp 
which is an odd function of m as required by time-reversal 
symmetry. The experimentally observed two consecutive 
metamagnetic transitions can be reproduced by a suit- 
able designed form of F{m) sketched in Fig. [T] a where 
only the part with m > is shown. It has two com- 
mon tangent lines marked with dotted and dashed lines, 
which touch the curve at points of mi. 2 and 7713^4, respec- 
tively. In the external magnetic field h, the solution of 
TO satisfies the equation: 



Therefore, the slopes of the two common tangent lines 
hi 2 can be interpreted as the fields at which magneti- 
zation jumps from mi to TO2 and from TO3 to m4, i.e., 
metamagnetic transitions occur. When h lies between 
these two transitions, the magnetization to evolves con- 
tinuously. 

The development of nematic orders between the two 
successive metamagnetic transitions is triggered by the 
(7(m)-term. Although the nematic instability is enhanced 
by Eq. [4l they are still a weaker instability compared 
to ferromagnetism (metamagnetism) . This is because 
the condition for Pomcranchuk instability in the d-wave 
channel in 2D, i.e., F2''^ < —2, is more stringent that in 
the s-wave channel, i.e. Fg < — 1. We assume that the 
charge and spin nematic channels are close to be critical 
but not yet, i.e., Vc^sp are small but positive. Due to the 
hybridization term of g{m), the eigen-order parameters 
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FIG. 2: (a) Dispersion of the 7+ band for k — (fc^,0) and 
{kx,n)- (b) The DOS of the 7± bands whose peaks corre- 
sponds to the van Hove singularities ai k — (0, 0), (0, tt), (tt, 0). 

arise from the diagonalization of the quadratic terms of 
^c^sp as. 

71+ — cos 07lc -f sin 07ls , 

rt_ = — sin 07ic + cos 07is , (8) 

where tanh20 = 2g{m)/{rc — 7-^^). The corresponding 
eigenvalues read: 

r± = \{rc^rsp±^ (r, - r^pY + 4[g(m)]2}. (9) 

The critical coupling for the 71+ channel to develop the 
instability is 

[,'(TO)]^ , kM]! > 1. (10) 

In the presence of the van Hove singularity of DOS, 
which is a common mechanism for meta-magnetism, all 
the parameters in the GL free energy could change signif- 
icantly so that the distribution of the dimensionless cou- 
pling function g' {m) may not be smooth. For the nematic 
order only occuring in the regime between two metam- 
agnetic transitions, 5' (to) must have a peak satisfying 
Eq. [TO] at TO2 < m < TO3 but is below the critical value 
elsewhere as sketched in Fig. [1] h. Roughly speaking, 
the underlying physics is that metamagnetism pushes the 
majority Fermi surface even closely to the van Hove sin- 
gularity which finally drives the nematic ordering^. The 
minority Fermi surfaces is pushed away from critical. In 
the following we will confirm this mechanism explicitly 
in a microscopic calculation. 

IV. MICROSCOPIC THEORY FOR THE 
BILAYER SrsRuiOj 

We next make the connection to Sr3Ru207 by exploit- 
ing a microscopic model with dxz and dy^-orbital bands. 
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They have the bonding and anti-bonding bands with a 
large bilayer splitting at the order of their band widths. 
Only the bonding bands, the even combination of the 
bilayer orbitals, are considered because their Fermi sur- 
faces are close to van Hove singularity which enhances 
interaction effects. Because of orbital hybridization, they 
form closed Fermi surfaces. As moving around this Fermi 
surface, the orbital configuration varies between d^z and 
dyz exhibiting a d-wave pattern, thus the mechanism il- 
lustrated in Sect. HIl applies. 

The band Hamiltonian Hq reads 



E 

ka 



|e td^ . d 



xz.ka 



krr yz,k<T 



h.c.)], 



(11) 



with 



^xzk ~ — 2i|| COS fcj, — 2tj^ COS /cj, — 4i' COS /cj; COS A:, 
e^^ ^: = — cos fca; — 2<ii cos /c„ — 4i' cos /c-K cos /c- 
Ar — 4t" cos fcr cos fc. 



1" 



V 



(12) 



i|| and t±_ are the nearest neighbor longitudinal and trans- 
verse hopping integrals for the d^z and dyz-orbitals; t' 
and t" are the next-nearest neighbor intra and inter- 
orbital hopping, respectively. The resulting diagonalized 
band Hamiltonian is: 



Ho 



E 

ka 



i- ' +ka '+kcr h I -ka 



ka 



(13) 



where 



k 2 



xz^k ' "yz 

and the band eigen-operators reads 

1+^Ua = cos</)^d^^^j^-Ksin</)^d^^j:^, 

^-ka ^ ^ 't^k'^xz ka + ^k'^vz ka 

with the hybridization angle satisfying 

—At" cos kx cos k. 



tan 2(/)r; 



2A£ 



yzM 



til 



t± COS kx — cos ky 



(15) 



(16) 



This band hybridization pattern is a lattice version of 
Eq. [2] with only the dx2-y2 channel but not the dxy chan- 
nel. We choose the parameter values of {t^^,t±,t' ,t") — 
(1.0, 0.145, 0.0, 0.3), and plot the density of states (DOS) 
in Fig [2l Two peaks in the DOS exist in both 7± 
bands, which correspond to the van Hove singularities 
at fe = (0,0),(0,7r),(7r,0). 

We take the general multi-band Hubbard model that 
are widely used in literatures for the interactions as 

Hint ^ U ^ na^{i)nai{i) + Vy^jixz{i)nyz{i) 

i,a—xz,yz i 



J 2_,{'^^2(*) ■ ^y^(^) ^ -^nxz{i)nyz{i)} 



+ ^Y.^dl^A^i)dl^^^{i)dyz^i{i)dyz,^{i) + h.c.}, 

(17) 

where Ua.a are particle number operators in orbital a with 
spin (t; Ua = ?T.a,T + "-a,!! are spin operators in orbital 
a. The [/-term is the intra-orbital repulsion; the y-term 
is the inter-orbital repulsion for the spin triplet configu- 
ration of two electrons; the J-term represents the Hund's 
rule physics; the A-term describes the inter-orbital pair- 
ing hopping. 

In the absence of the orbital hybridization, it is not 
conclusive that the Hubbard model can give rise to ne- 
matic transitions because on-site interactions only con- 
tribute to the s-wave channel. We define the charge 
nematic order 2nc — n^z — nyz, the spin nematic or- 
der Ugp — Sxz — Sy^, and the ferromagnetic order m — 
Sxz + Sy^. The mean- field theory in the eigen-basis of 7± 
reads: 



(18) 



^ ^ ' aka 

ka.a—zL 



where 



U J 



U J 



J U 



yrn = ^ + -,V,,^^--,V,^V+-~ -, (19) 



and the less important inter-band coupling terms 

'+ka 

cle spectra become: 



7]^- 7_feCT dropped. The renormalized single parti- 



^aa = ES - ^-aVmm-aiycnc + crVspnsp) cos24>k, (20) 

where (t)k are the hybridization angle defined in Eq. 1161 
The self-consistent equations for the order parameters 
read: 

= -^^{n/,+,^(fc) -n/,_,<,(fc)}cos20g. 



k.a 



nsp = — ^a{n/,+,o.(fc) -n/_,<,(fc)}cos20g, 

k.a 



k,i7 

where the cos2(/)j factor represents the dx2-y2 symmetry 
of the nematic orders of ric and rigp. 

We are ready to study the zero-temperature phase di- 
agram with the Zeeman energy term as: 



a—xz^yz 



(22) 



/i is tuned to reach the filling of n = 3.48 at /i = in 
the bilayer bonding bands of dxz and dyz so that the FS 
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FIG. 3: m and the nematic order parameter for majority 
spin band (nc| = (nc + nsp)/2) versus the external field h. 
The topologies of the FS of the 7+ band in each phases are 
sketched in the insets, where the solid and dashed lines rep- 
resent the FS of the majority and the minority spin bands. 

of the 7+ band is very close to the van Hove singularity 
while that of the 7_ band is far away from the singularity. 
The interaction parameters are chosen as {U, V, J, A) = 
(2.7, 2.5, 0, 0), which remarkably reproduce the isotropic- 
nematic-isotropic phase transition as a function of h as 
we expect from the Ginzburg-Landau analysis. 

The magnetization m and the nematic order of the ma- 
jority spin Fermi surface nc,T = |("-c + ngp) are depicted 
in Fig. [3l while that of the minority spin Fermi surface 
is much smaller (not shown) because its Fermi surface is 
pushed from the van Hove singularity. As explained in 
Ref.— , the first jump of m distorts the Fermi surface to 
touch the van Hove singularity points along one of the 
ab axes, and the second jump restores the four fold sym- 
metry of Fermi surface to cover the singularity in both 
directions. Compared to Refi^, the range of h for the 
nematic state is significantly reduced in agreement with 
experiment due to inclusion of the magnetic order m self- 
consistently in the solution. Furthermore, the first jump 
of m is larger than the second one as consistent with 
the experiments. This feature can be understood by the 
asymmetric profile of the DOS at the van Hove singular- 
ity which drops faster in the higher energy side as shown 
in Fig. O After the first transition, part of the FS of 
the majority spin has moved beyond the van Hove sin- 
gularity. This reduces the DOS at the second jump and 
leads to a weaker second jump where the nematic order 
disappears. 

We learned another independent and beautiful work 
by Raghu et al.'^^ which was posted at the same time as 
ours. The same mechanism is proposed for the nematic 



state observed in Sr2Ru307 based on the quasi-one di- 
mensional bands. They used a more detailed band struc- 
ture and further considered the spin-orbit coupling effect. 
We further performed the Ginzburg-Landau analysis for 
the competition between ferromagnetization and nematic 
ordering. 



V. CONCLUSION AND OUTLOOK 

In this article, we have studied the Fermi liquid in- 
stability of the Pomeranchuk type in orbital band sys- 
tems. Orbital band hybridization significantly enhances 
the Landau interaction functions in high partial wave 
channels, thus providing a mechanism for the nematic 
states or unconventional magnetism from conventional 
interactions. Consequentially, metamagnetism (ferro- 
magnetism) induces the nematic behavior even with the 
onsite multi-band Hubbard interactions. 

This mechanism is applied to the t2g system of 
Sr2Ru307 by attributing the observed nematic behav- 
ior to the hybridized quasi-one dimensional bands of d^z 
and dyz, which is the major difference between our work 
and Refiii^i^. Many open questions still need future ex- 
ploration. In particular, the quick suppression of the ne- 
matic behavior with the in-plane magnetic field might 
result from the orbital effect due to the bilayer splitting 
as pointed out in Ref.^^, or from the spin-orbit coupling 
effect. 

The mechanism presented in this article also very gen- 
eral. It is essentially a Berry phase effect which naturally 
arises from electron systems with non-trivial band struc- 
tures, such as spin-orbit coupling system and graphene. 
It can be understood as a conventional interaction ac- 
quires a non-trivial nature after projected onto a non- 
trivial band structure. We further predict that the ne- 
matic ordering arises at the ferromagnetic transition in 
spin-orbit coupling systems as a result of the hybridiza- 
tion between two spin components in the band struc- 
ture. For example, the in-plane ferromagnetic ordering 
in the quasi-2D Rashba systems and the surface states 
of the topological insulators induces the p-wave distor- 
tion of Fermi surfaces, which will be presented in a later 
publication^. 
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